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High  fidelity  simulation  of  atomization  in  diesel 

engine  sprays 

By  L.  Bravof,  C.  B.  Ivey,  D.  Kim|,  and  S.  T.  BoseJ 


A  high  fidelity  numerical  simulation  of  jet  breakup  and  spray  formation  from  a  complex 
diesel  fuel  injector  has  been  performed.  A  full  understanding  of  the  primary  atomization 
process  of  diesel  fuel  injection  has  not  been  achieved  for  several  reasons,  including  the 
difficulties  in  accessing  the  optically  dense  region.  Due  to  recent  advances  in  numerical 
methods  and  computing  resources,  high  fidelity  simulations  of  realistic  atomizing  flows 
are  currently  feasible,  providing  a  new  mechanism  to  study  the  jet  breakdown  process. 
In  the  present  study,  a  novel  volume-of- fluid  (VOF)  method  coupled  to  a  stochastic  La- 
grangian  spray  (LSP)  model  is  employed  to  simulate  the  atomization  process.  A  common 
rail  fuel  injector  is  modeled  by  a  nozzle  geometry  provided  by  the  engine  combustion  net¬ 
work  (ECN).  The  working  conditions  correspond  to  a  single  90  /im  orifice  JP-8  fueled 
injector  operating  at  90  bar  and  373  K  and  releasing  into  a  100%  nitrogen,  29  bar,  300 
K  ambient  with  a  Rei  =  16,071  and  Wei  =  75,334,  putting  the  spray  in  the  full  at¬ 
omization  mode.  The  experimental  dataset  from  Army  Research  Lab  (ARL)  is  used  for 
validation  and  the  Kelvin-Helmholtz/Rayleigh-Taylor  (KH-RT)  breakup  model  (Reitz  & 
Bracco  1979)  is  used  for  verification,  both  in  terms  of  spray  angle.  Droplet  distributions 
of  the  simulated  spray  are  provided  for  future  experimental  comparisons  and  secondary 
atomization  simulations  using  LSP  modeling. 


1.  Introduction 

To  date,  one  of  the  main  bottlenecks  in  engineering  spray  modeling  of  combustion 
systems  is  an  accurate  description  of  the  primary  atomization  process.  Several  contem¬ 
porary  numerical  solvers  adopt  coarse  approximations  in  the  dense  region  based  on  in¬ 
jecting  nozzle-sized  particles  subject  to  Kelvin-Helmholtz  type  instabilities  coupled  with 
Lagrangian  particle  tracking  techniques  (Som  &  Aggarwal  2009;  Senecal  2012).  Signif¬ 
icant  success  has  been  achieved  with  these  methods;  however,  they  require  a  posteriori 
knowledge  of  the  spray  process  to  calibrate  the  model,  making  an  a  priori  investiga¬ 
tion  impossible.  Historically,  the  advancement  of  primary  atomization  models  has  been 
hindered  by  the  well-known  difficulties  in  measuring  the  optically  dense  spray  region.  Al¬ 
though  experimentalists  have  had  success  with  modern  methods,  such  as  ballistic  imaging 
and  x-ray  techniques,  extraction  of  full  four  dimensional  information  with  sufficient  spa¬ 
tial  and  temporal  resolution  for  a  detailed  analysis  is  still  infeasible  (Linne  et  al.  2006; 
Wang  et  al.  2006;  Coletti  et  al.  2014). 

Remarkable  progress  has  been  made  in  recent  years  in  the  development  of  robust  nu¬ 
merical  methods  for  handling  interfaces,  enabling  researchers  to  perform  highly  resolved 
simulations  of  multiphase  flow  (Gorokhovski  &  Herrmann  2008).  Desjardins  reported 
on  the  development  of  the  level  set /ghost  fluid  method  utilizing  high  order  schemes  to 
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study  multiphase  flows  (Desjardins  et  al.  2008).  The  accuracy  of  the  numerical  technique 
was  corroborated  through  several  studies,  including  the  atomization  of  a  liquid  diesel 
jet  at  moderately  low  Reynolds  number,  Re  =  3000.  In  a  related  study,  Desjardins  con¬ 
ducted  detailed  numerical  simulations  (DNS)  of  primary  atomization  for  several  values  of 
Reynolds  and  Weber  numbers  (2000  <  Re  <  3000,  500  <  We  <  2000),  reporting  on  ve¬ 
locity  statistics  across  the  turbulent  jet  (Desjardins  &  Pitsch  2010).  More  recently,  using 
a  refined  level  set  grid  approach,  Herrmann  discussed  the  impact  of  finite  grid  resolu¬ 
tion  on  the  phase  interface  geometry  of  the  liquid  jet  core  under  diesel  engine  conditions 
with  Re  =  5000  and  an  injection  velocity  of  100  m/s  (Herrmann  2011).  In  this  work,  it 
was  reported  that  turbulence  is  the  driving  mechanism  of  atomization  within  the  first 
20  diameters  downstream  of  the  injector.  It  was  also  determined  that  6  grid  points  are 
needed  to  obtain  grid-independence  of  larger  scale  drops.  These  studies  provide  a  criti¬ 
cal  database  to  drive  the  next-generation  spray  model  development.  Note  that,  as  with 
most  other  DNS  studies,  no  quantitative  comparison  to  experimental  data  is  typically 
provided. 

The  need  to  accurately  model  two  phase  atomizing  flows  in  high  speed  jets  is  partic¬ 
ularly  important  in  diesel  injectors  where  the  quality  of  the  fuel  and  oxidizer  mixing  is 
essential  for  lean  combustion.  Fuel/air  mixture  formation  is  also  a  very  important  factor 
in  increasing  engine  efficiencies  and  power  densities.  Spray  and  atomization  characteris¬ 
tics  have  to  be  considered  to  optimize  the  design  of  the  combustion  chamber  to  reduce 
exhaust  emissions  and  to  improve  combustion  performance.  Also,  diesel  spray  charac¬ 
teristics  can  be  influenced  by  the  injector  geometry,  the  injection  parameters,  and  the 
flow  mixing  inside  the  combustion  chamber.  Therefore,  simulations  should  account  for 
system  level  complexities,  including  real  injector  features  to  accurately  predict  realistic 
spray  dynamics. 

Additionally,  the  computational  expense  of  resolving  all  the  critical  length  scales  at 
large  Weber  number  is  prohibitively  large,  so  the  application  of  DNS  to  simulate  realistic 
diesel  injector  conditions  has  been  severely  limited.  A  liquid  jet  moving  with  an  0(100) 
m/s  relative  velocity  with  respect  to  the  quiescent  gas  can  generate  droplets  with  diam¬ 
eters  as  a  small  as  a  few  microns  (Desjardins  &  Pitsch  2010).  Hence,  there  is  a  need  to 
develop  more  accurate  engineering  breakup  spray  models  for  the  primary  and  secondary 
breakup  modes  to  reduce  the  computational  cost  when  simulating  a  spray  filled  domain. 

The  objective  of  the  current  work  is  to  investigate  the  atomization  behavior  of  a 
high-speed  single  hole  jet  with  complex  internal  geometry.  A  novel  unstructured  volume- 
of-fluid  (VOF)  method  has  been  adopted,  which  is  geometric  and  un-split,  enforcing 
exact  mass  conservation  on  an  unstructured  grid  (Kim  et  al.  2013,  2014).  In  this  study, 
the  VOF  method  is  coupled  to  the  Lagrangian  spray  (LSP)  framework  to  increase  the 
computational  efficiency  and  to  apply  sub-grid  atomization  models.  Hence,  the  phase 
interface  resolved  by  the  grid  is  captured  by  the  VOF  method,  while  the  under  resolved 
small  scale  droplets  are  transferred  from  the  VOF  interface  representation  to  the  LSP 
particle  tracking  and  further  breakup  is  handled  by  a  stochastic  breakup  model.  Mea¬ 
surements  were  conducted  at  the  Spray  and  Combustion  Research  Facility  of  the  Army 
Research  Lab  (ARL)  to  qualitatively  complement  the  simulation  and  to  validate  the  sim¬ 
ulation  in  terms  spray  angle.  The  Kelvin-Helmholtz/Rayleigh- Taylor  (KH-RT)  breakup 
model  (Reitz  &  Bracco  1979)  is  used  to  verify  the  simulation  in  terms  of  spray  angle. 
Droplet  distributions  were  generated  for  future  experimental  comparisons  and  secondary 
atomization  models. 
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2.  Methods 

2.1.  Simulation 

The  detailed  numerical  simulation  of  the  interface  was  performed  using  a  novel  geometric 
unsplit  VOF  method  that  is  conservative  on  unstructured  meshes  coupled  to  a  stochastic 
Lagrangian  spray  (LSP)  framework.  The  geometric  VOF  method  ensures  discrete  con¬ 
servation  and  boundedness  of  the  volume  fraction,  F,  by  utilizing  non-overlapping  flux 
polyhedra  for  donor  volumes  (see  Ivey  &  Moin  2012  for  a  description  of  flux  polyhedra). 
The  volume-of-fluid  method  uses  piecewise  linear  interface  calculation  (PLIC)  represen¬ 
tation  of  the  interface,  requiring  an  interface  normal,  n.  h  is  calculated  from  an  auxiliary 
level  set,  G,  that  was  updated  using  the  geometric  advection  algorithm  to  keep  it  consis¬ 
tent  with  F.  Curvature,  k  is  also  calculated  from  the  G  using  the  direct  front  curvature 
method  (Herrmann  2006).  After  n  and  k  are  calculated,  G  is  reconstructed  to  strictly 
follow  the  PLIC  representation  using  a  bisection  algorithm  to  calculate  the  local  G  field 
(to  enforce  F)  and  the  reconstructed  distance  function  (Cummins  et  al.  2005)  method 
to  propagate  the  G  field  throughout  the  band.  For  consistency  (and  stability),  mass  and 
momentum  are  convected  using  the  geometric  VOF  method.  To  diminish  the  strict  over¬ 
flow  time  step  requirements  on  VOF  schemes,  multiple  frozen  velocity  advection  updates 
(taken  as  3  here)  are  performed  for  each  momentum  step.  Several  validation  studies  have 
been  performed  that  tested  the  numerical  accuracy  and  robustness  of  the  solver  (Kim 
et  al.  2013,  2014). 

For  computational  efficiency,  the  VOF  representation  is  coupled  to  a  stochastic  La¬ 
grangian  spray  (LSP)  representation.  Interfacial  features  are  characterized  by  contiguous 
regions  of  F  >  0.  Under  resolved  (taken  <  53  cells  here)  interfacial  flow  structures  are  re¬ 
placed  by  a  spherical  droplet  of  equal  volume.  No  additional  breakup  model  is  employed 
here  so  as  to  capture  the  instabilities  leading  to  atomization  directly.  LSP  droplets  fol¬ 
low  particle  drag  laws  and  a  stochastic  breakup  model  based  on  critical  Webber  number 
(Ham  et  al.  2003). 

The  simulation  was  conducted  on  the  nozzle  geometry  available  through  the  engineer 
combustion  network  (ECN)  with  flow  conditions  corresponding  to  an  injection  pressure 
of  90  bar,  a  background  pressure  of  20  bar,  and  a  bulk  jet  exit  velocity  of  127  m/s.  The 
Reynolds  and  Weber  number  were  calculated  with  JP-8  fuel  properties  database  (at  373 
K)  and  yielded  a  value  of  Rei  =  16,071  and  Wei  =  75,334.  The  critical  length  scales, 
determined  from  the  problem  configuration,  range  from  the  nozzle  orifice  integral  scale 
(Zj  =  d  =  90  pea),  to  the  viscous  scale  (0.09  pin),  and  down  to  Kolmogorov  critical  radius 
(0.2  pm).  The  pipe  flow  viscous  scale  and  the  critical  radius  are  defined  as, 

lv  =  1/0.2  Re"7/8  and  lcr  =  .  (2.1) 

A  large  eddy  simulation  dynamic  Smagorinsky  model  was  adopted  to  treat  the  smallest 
flow  structures  and  for  computational  efficiency.  In  the  nozzle,  a  wall  resolved  approach 
was  utilized  with  Ax+,  A y+  ~  1  near  the  wall  (calculated  from  lv),  Ax+,  A y+  ~  50  near 
the  center,  and  A z+  ~  70  (see  Figure  1).  The  mesh  refinement  is  cartesian  in  nature, 
forcing  Ax  =  Ay  =  A r.  A  mesh  of  ~  60  million  hexahedral  control  volumes  was  used  to 
discretize  the  domain.  Mesh  points  were  concentrated  in  the  jet  spray  envelope  to  resolve 
the  interface,  while  coarse  buffer  regions  were  added  in  the  radial  and  axial  regions  in 
order  to  reduce  the  impact  of  the  boundary  conditions  (see  Figure  2).  Within  the  spray 
envelope,  A/lcr  ~  0.5  near  the  unperturbed  nominal  interface  to  capture  the  instabilities 
and  A /lcr  ~  40  throughout  the  rest  of  the  spray  envelope  to  leverage  the  LSP  framework 
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for  savings.  The  numerical  boundary  conditions  in  the  ambient  region  were  a  constant 
inflow  set  to  5%  the  bulk  jet  velocity,  a  slip  condition  over  the  radial  boundaries,  and  a 
uniform  outflow. 


2.2.  Experiment 

The  experiments  in  this  study  were  conducted  at  the  Army  Research  Laboratory  by 
injecting  a  high  speed  JP-8  fueled  spray  into  a  high  temperature  pressure  (HTPV)  flow 
through  chamber.  The  HTPV  is  designed  to  reach  a  maximum  pressure  of  150  bar  and 
a  maximum  temperature  of  1000  K  using  a  BOSCH  common  rail  injection  system  for 
precise  fuel  delivery  (Kurmann  et  al.  2014).  The  vessel  is  equipped  with  closed  loop 
control  for  pressure  and  temperature.  The  flow  through  chamber  is  held  constant  at  58 
m3/hr.  An  on-site  nitrogen  generator  produced  the  necessary  nitrogen  for  testing,  which 
was  maintained  at  99%  purity  during  experiments.  To  allow  for  optical  access,  the  vessel 
is  outfitted  with  3  fused  silica  windows  with  dimensions  of  147  mm  diameter  by  85  mm 
thickness.  To  protect  the  85  mm  thick  pressure  windows  from  fuel  contamination,  6  mm 
thick  fused  silica  windows  are  placed  between  the  85  mm  windows  and  the  spray  zone. 

High  speed  near  field  spray  region  images  were  acquired  using  a  single  LED  light  source 
and  a  Photron  SA5  camera  operating  at  90,000  fps  for  line  of  sight  measurements.  For 
the  experiments  presented,  image  size  was  set  to  320  by  192  pixels  and  the  corresponding 
scaling  was  5.6  ^m/pixel.  Chamber  conditions  were  set  to  20  bar  and  300K  prescribing  a 
density  ratio  of  34  kg/rn3.  Fuel  injection  pressure  was  set  to  90  bar  with  a  3  ms  injection 
duration,  and  a  total  of  2.2  mg  of  injected  mass  was  measured  via  an  IAV  injection 
analyzer.  Figures  3  and  4  and  show  two  instances  where  the  spray  behavior  goes  from 
transitional  to  fully  atomizing  mode. 


2.3.  Theory 

Reitz  dispersion  model  has  been  proposed  to  study  the  spray  angle  by  employing  aero¬ 
dynamic  arguments  (Reitz  &  Bracco  1979).  It  includes  the  ratio  of  the  Reynolds  and 
Weber  number  of  the  liquid  flow  in  the  function  /(y)  and  it  is  written  as 

=  I  fe) /h)- 


(2.2) 


Figure  3.  High  speed  near  field  image  of 
JP-8  spray  during  transition  state 


Figure  4.  High  speed  near  field  image  of 
JP-8  spray  during  fully  atomized  state 


Figure  5.  Q  criterion  isosurface  colored  by  Figure  6.  U  contour  slice  at  jet  centerline, 
streamwise  velocity  in  the  diesel  spray  injector 
as  viewed  from  the  nozzle  exit. 


where  pg  and  pi  are  the  liquid  and  gas  density,  and  A  is  a  constant  that  depends  on  the 
nozzle  design,  A  =  3.0  +  0.28Zo/do;  with  d0  is  the  nozzle  diameter  and  l0  is  the  length  of 
the  nozzle  hole.  The  parameter  /(y)  is  function  of  the  physical  properties  of  the  liquid  and 
injection  velocity  and  defined  as,  /(y)  =  V3/6[l  — exp(— lOy)]  with  7  =  (Rei/Wei)2 pi/ pg 
where  the  injection  velocity  is  based  on  Bernoulli  arguments,  Vinj  =  Cv\j2tAP / p. 


3.  Results 

The  nozzle  flow  turbulence  was  visualized  by  sampling  the  velocity  flow-field  and  using 
the  classical  Q  criterion  defined  as  Q  =  l/2(||fijj||  —  ||<SV/||)  and  colored  by  the  streamwise 
velocity  component,  U  (see  Figure  5).  The  Q  criterion  isosurface  shows  hairpin  vortex 
structures  arising  from  the  interaction  of  the  fluid  with  the  wall,  having  peak  streamwise 
velocity  magnitudes  towards  the  center  of  the  pipe  flow.  The  vortex  structures  are  ir¬ 
regular  as  compared  to  a  traditional  pipe  flow;  this  discrepancy  can  by  explained  by  the 
favorable  pressure  gradient  of  the  nozzle  and  the  lack  of  perfect  symmetry  in  the  exper¬ 
imentally  modeled  geometry.  Nevertheless,  the  emerging  flow  field  can  be  characterized 
as  turbulent  and  the  resultant  jet  breakdown  can  be  interpreted  as  being  in  the  spray 
atomization  regime.  The  turbulent  inflow  combined  with  the  jet  instabilities  lead  to  the 
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Figure  7.  F  isosurface  and  Lagrangian  particles  alongside  the  diesel  spray  injector. 


chaotic  jet  behavior  illustrated  by  the  streamwise  velocity  contour  in  Figure  6. 

Figure  7  shows  the  spray  formation  process  at  steady  full  atomization  conditions  as 
visualized  using  the  F  =  0.5  isosurface  and  the  LSP  tracked  droplets  with  reference 
to  the  nozzle  geometry.  Note  the  growth  of  the  hydrodynamics  instabilities  and  the 
resultant  spray  cone.  The  present  simulation  solely  models  the  static  fully  opened  valve 
configuration. 

Figures  8-11  show  the  radial  variations  of  U  and  F  statistics  with  axial  distance.  Fig¬ 
ure  8  shows  the  use  of  5%  co-flow  field  used  to  stabilize  the  solution  (U  ^  OVr  >  D/ 2 
)  and  the  experimental  bulk  velocity  (at  0  D).  Indicative  of  growing  jets,  the  U  pro¬ 
files  broaden  downstream.  Figure  9  shows  the  mean  F  distribution,  which  includes  the 
equivalent  volume  fraction  of  the  Lagrangian  particles.  As  in  [7,  but  to  a  smaller  de¬ 
gree,  F  profiles  broaden  downstream.  The  F  profiles  decrease  in  height  downstream 
demonstrating  the  breakdown  and  fluctuation  of  the  liquid  jet.  Note  that  the  dispersion 
characteristics  between  velocity  and  volume  fraction  fields  are  quite  different,  demon¬ 
strating  the  entrainment  effect  on  the  velocity  fields  and  the  conservation  of  mass  of  the 
volume  fraction  fields.  The  turbulent  kinetic  energy  (TKE)  radial  profiles  are  shown  in 
Figure  10.  The  nozzle  pipe  flow  injects  sharp  TKE  peaks,  generating  turbulent  structures 
at  fluid  interface.  The  TKE  facilitates  the  jet  breakdown  and  grows  with  the  shear  layer. 
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Figure  8.  Mean  streamwise  velocity  as  a 
function  of  radius  for  5  downstream  stations: 
z/D  =  { 0,4,8,12,16}  (dark  to  light). 


Figure  9.  Mean  volume  fraction  (including 
Lagrangian  spray)  as  a  function  of  radius  for 
5  downstream  stations:  z/D  =  {0,4,8,12,16} 
(dark  to  light). 


r /D 

Figure  10.  Mean  turbulent  kinetic  energy  as 
a  function  of  radius  for  5  downstream  stations: 
z/D  =  {0,  4, 8, 12, 16}  (dark  to  light). 
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Figure  11.  Root  mean  square  of  volume  frac¬ 
tion  (including  Lagrangian  spray)  as  a  func¬ 
tion  of  radius  for  5  downstream  stations: 
z/D  =  { 0,4,8,12,16}  (dark  to  light). 


The  intensity,  F  root  mean  square  (RMS)  profiles,  in  Figure  11  shows  the  intensity  fields 
increasing  with  the  growth  of  the  shear  growth  layer. 

The  statistics  have  not  fully  converged,  as  shown  in  the  Figure  8,  so  in  order  to  extract 
a  preliminary  estimate  of  the  spray  angle,  0,  we  fit  a  Gaussian  curve  to  each  of  the  mean 
velocity  profiles.  Figure  12  shows  the  Gaussian  fits,  which  are  not  centered  at  the  origin, 
evidencing  the  lack  of  axisymmetry  of  the  real  geometry  (the  spray  comes  off  at  a  angle) . 
The  full  width  half  maximum  (FWHM)  of  the  Gaussian  fits  to  the  mean  streamwise 
profiles  is  shown  in  Figure  13;  the  change  of  the  linear  fit  to  the  FWHM  with  axial 
distance  provides  the  spray  angle,  8. 

The  experimental  spray  angle  was  determined  via  line  of  sight  observations  of  200 
spray  images  (sampling  frequency  of  11.1  ps)  while  tracking  the  interface  of  the  jet  core 
region  with  respect  to  the  jet  centerline.  The  total  sampling  time  corresponds  to  2.2  ms 
in  the  full  atomization  spray  mode  (Re  =  16,  071  and  Oh  =  0.017).  The  initial  transients 
were  not  included  in  the  procedure  for  consistency  with  the  simulation  results.  The  spray 
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Figure  12.  Gaussian  fit  to  mean  streamwise 
velocity  as  a  function  of  radius  for  5  down¬ 
stream  stations:  z/D  =  {0,4,8,12,16}  (dark 
to  light). 


Figure  13.  Full  width  half  maximum  of  gaus- 
sian  fit  to  mean  stream  wise  velocity  as  a  func¬ 
tion  of  downstream  location.  8  is  the  angle 
of  linear  fit  to  the  full  width  half  maximum 
(shown  in  gray). 


angles  extracted  from  Reitz  theory,  simulation,  and  experiment  are  respectively 

0Reitz  =  5.8325°  ,  0sim  =  5.1125°  and  9exp  ~  4°.  (3.1) 

The  simulation  and  theoretical  dispersion  results  appear  to  be  in  good  agreement. 
However,  there  is  a  1°  discrepancy  with  the  experiments  featuring  a  spray  angle  of  4°. 
The  differences  can  be  in  part  due  to  variations  in  the  nozzle  orifice  diameter  and  nozzle 
shape  that  arise  from  manufacturer  fabrication  eccentricities.  This  has  been  reported  in 
the  literature  previously  for  ECN  type  injectors  where  geometric  inconsistencies  were 
thoroughly  reported  (Kastengran  et  al.  2012).  The  impact  of  these  discrepancies  can 
clearly  affect  the  spray  parameters. 

Secondary  atomization  models  require  an  initial  droplet  spray  profile,  parameterized 
by  droplet  diameter  and  distance  away  from  jet.  Droplet  sizes  and  counts  are  calculated 
from  the  combination  of  the  LSP  particles  and  the  resolved  VOF  features  (calculated 
using  the  same  method  to  transfer  from  VOF  to  LSP  for  under  resolved  features,  but 
with  a  larger  domain  cell  count).  Figure  14  shows  the  average  number  of  droplets  for  a 
given  droplet  diameter.  The  droplets  range  from  approximately  1  —  10  /im,  the  average 
droplet  diameter  is  2.9786  /im,  and  the  most  likely  droplet  diameter  is  1.5055  /im.  Droplet 
diameters  <  1  /im  were  spontaneously  evaporated  and  were  not  tracked.  Figure  15  shows 
the  average  droplet  diameter  at  a  given  distance  away  from  the  jet  center.  The  average 
droplet  size  increases  away  from  the  jet  going  from  approximately  2  —  7  /im. 


4.  Conclusions 

In  this  investigation  a  high  fidelity  simulation  approach  was  adopted  to  study  the  atom¬ 
ization  physics  of  a  diesel  injector  with  detailed  nozzle  internal  geometry.  The  nozzle  flow 
field  was  characterized  through  visualizations  of  Q  isosurfaces  for  turbulence  patterns. 
The  complexity  of  the  geometry  and  system  dynamics  was  characterized  by  a  snapshot  of 
the  volume  fraction  isosurface  and  Lagrangian  droplets.  Also,  mean  streamwise  velocity 
and  volume  fraction  statistics  show  the  structure  of  the  high  speed  jet.  The  turbulent 
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Figure  14.  Number  of  droplets  as  a  function 
of  diameter 


Figure  15.  Mean  droplet  diameter  as  a 
function  of  radius 


kinetic  energy  and  volume  fraction  inensitity  profiles  characterize  the  interfacial  mixing 
processes.  Comparison  with  Reitz  spray  theory  and  ARL  measurements  of  the  near  noz¬ 
zle  flow  field  show  that  the  simulation  captures  the  correct  dispersion  characteristics. 
The  spray  was  further  characterized  using  droplet  size  and  spatial  distribution  plots. 
Further  work  is  presently  on  the  way,  using  higher  resolution  to  establish  numerical  con¬ 
vergence  and  to  capture  the  hydrodynamic  flow  instabilities  for  comparison  with  classical 
instability  models. 
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